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ABSTRACT 


This  r#p*rtis  concerned  with  the  approximation  of  acoustical  noise 
fields  by  the  use  of  modulated  pure  tone  sources.  Given  a  specific  random 
signal  the  methods  of  orthogonal  expansion  and  Fourier  expansion  are  studied. 
The  integral  equation  arising  from  the  orthogonal  expansion  is  solved  for 
the  normal  functions  when:  a)  the  power  spectrum  is  given  as  the  ratio  of 
two  polynomials,  b)  Markoff  Autocorrelation  function,  c)  stationary  band 
limited  Gaussian  white  noise  of  mean  zero,  d)  band  limited  stationary  normal 
white  noise  on  a  finite  interval,  e)  white  noise.  Deterministic  and  random 
frequency  and  amplitude  modulation  and  random  switching  are  examined  as 
ways  to  broaden  spectral  components. 
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I.  Introduction 


Because  of  the  many  applications,  it  is  found  useful 
to  represent  a  random  process  in  a  suitable  series  of 
particular  functions  over  some  finite  or  infinite  interval. 
Considerable  study  has  been  devoted  to  the  problem  of 
expansions  (Ref.  1,2,3)*  Orthogonal  expansions  have 
received  particular  attention  (Ref.  2). 

The  present  report  deals  with  some  aspects  of  the 
expansion  of  random  functions  into  orthogonal  series. 
Consideration  is  also  given  to  how  to  broaden  the  spectrum 
of  a  single-tone  signal  by  modulation  of  various  sorts. 

Some  attention  is  also  given  to  optimal  approximations  of 
various  statistical  characteristics  of  actual  random  noise 
fields. 


Manuscript  released  by  author  October  1,  1964,  for  publication  as 
an  DTD  Technical  Documentary  Report. 
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II •  Representation  of  Acoustic  Noise  Fields 

a)  Orthogonal  Expansion  of  Random  Functions 

It  is  desired  to  represent  a  randomly  varying  sound 
pressure  p(t)  (see  Figure  1)  as  a  series  of  orthogonal 
functions  with  orthogonal  random  variables  as  coefficients. 


P(t) 


Figure  1#  Topical  Recording  of  Random  Sound  Field 
In  particular,  we  require 


(1) 


(2) 


To  achieve  these  results,  consider 


P(t)  p*(f)tn(f)  -Xlx.l  1^1 


2 


Integrating  both  sides  of  this  result  from  t-  0  to  T  and  using 
the  second  of  (2),  ve  find 

j^*p(t)P*(f)tn(f)df .  v  ix.i  i\iv.(t>  Vk 

•1  l\l  'XA(t)  V» 

s 

tricing  expection  and  using  the  first  of  (2)  gives 

/ *  R(t,t'W  (t')dt'  -  |X  |2  t  (t) 
n  n  n 

o<  t  <  T  (5) 

R(t,f)  -  I(p(t)  p(t')J 

Thus  ve  can  have  such  an  expansion  as  (2)  provided  tQ(t)  are 
solutions  of  (3)  and  X  is  the  eigenvalue  associated  vlth 

Prom  (1)  and  (2)  it  follows  that 

-  irri\<t)  »<*>«  <x> 

Therefore  if  we  knew  the  eigenfunctions  t^(t),  and  the 
eigenvalues  X  ,  and  possessed  a  realized  value  of  p(t), 
ve  could  calculate  the  realized  values  of  the  orthogonal 
random  variables  X  •  All  that  is  needed  to  find  X  and 
f  (t)  is  the  autocorrelation  function  R(t,t' ).  This  is 
either  given  or  is  determined  from  an  actual  recording 
of  p(t). 

If  p(t)  is  a  gauss lan  random  process,  the  random 
variables  X  will  not  only  be  orthogonal  but  will  be 
independent  random  variables  with  a  normal  probability 
distribution.  The  problem  of  determining  the  probability 
distribution  function  for  p(t)  when  those  of  X^  are  given 
is  an  extremely  difficult  one  and  has  been  solved  only  in 
the  case  of  gauss  lan  variates.  It  is  to  be  noted  from  (4) 
that  if  p(t)  is  an  Independent  process,  X  will  be  a 
normal  variate.  One  can  easily  prove  this  by  application 
of  the  central  limit  theorem  of  probability  theory. 
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b)  Fourier  Expansion  of  BaaOem  Function* 

Because  of  the  difficulty  of  solving  the  integral  equation 
(3)  and  because  of  simplicity,  random  functions  are  expanded 
into  Fourier  series  in  most  practical  situations.  Unfortunately, 
the  random  coefficients  are  not  orthogonal  on  (0,T)  except  for 
periodic  processes.  We  can  show,  however,  that  the  coefficients 
are  asymptotically  orthogonal  as  the  length  of  the  recording 
increases . 

Consider  a  stationary  random  function  p(t)  defined  on  the 
interval  (0,T) .  The  Fourier  expansion  of  this  process  on  (0,T) 


is 


p<*)  mL  co*  V  +  bk  B±a  Vs1,  (5) 

k-0 


where 


T 


cos  k»  t  cos  t'dtdt' } 


o  o 


cos  ku  t  CM  iw  t'dtdt' 
o  o 


T 


Changing  variables  the  double  integral  beccnes 


cos  kw  tdt 


o 


Let  A  «  t/T,  then 


k 


(\\>  -  { 


kl  cos  fcnkXdX  f  r(t)  cos  2rJ(-  ♦  X)  dr 


T(l-X) 


At 


r1  r*0--*)  27rlT 

»  t  /  cos  2flkXdX  (cos  2rriX  /  R(t)  cos  — —  dr 

r*1-*)  2wJtx 
-  sin  2 ttJX  /  R(t)  sin  — —  dr) 


As  T  +  m,  X  j  0 


IE  (a^a^}  J'  cos  2i7kXdX  (cos  2iriX*  ~  S(*s^ )  + 


♦  sin  2iriX;  (0)  } 


(7) 


2tJ 


S(w^)  ■  Spectral  density  of  p(t)  and  «  - 

unless  i  is  chosen  such  that  I  In  w  >  0. 

i,T  -»  •• 

(7)  nay  he  evaluated  to  give 


u,  ->  0  aa  T  -*  «• 


“(\V  -.3^) 

Hie  normalised  ■  a ^<*^>*  *j  “  *1 

obey 


(8) 


E(Vi}  •*  >  k*'> 1 


(9) 


In  like  wanner  one  can  shoe  that 


«\V 


E(S  S  )  -♦  8 

Yi 1  ki 


,(8oV-*1 

E(loV  ‘,8ol 


k,i  >  1 


(10) 


Thus  the  expansion  coefficients  are  asymptotically  orthogonal. 
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It  is  of  value  to  compare  various  Fourier  expansions  for 
different  lengths  of  T  of  recordings  of  p(t)  vith  the  exact 
orthogonal  expansion  on  (0,T).  This  would  give  some  idea  of 
the  length  of  record  to  he  taken  so  that  a  Fourier  expansion 
would  he  as  good  as  the  exact  one. 

HI.  Solution  of  the  Integral  Equation  Arising  in  tig  Orthog^i 
Expansion  of  a  Random  Function. 

a)  General  Results  for  Processes  with  Rational  Spectral 
Density  Functions 

We  confine  ourselves  to  stationary  random  sound  fields  whose 
spectral  density  functions  are  rational  functions  of  frequency. 

In  this  case  the  Integral  equation  becomes 

JQ  *(*-*')  tn(t')dt’  -  l\|2*n(t),  0  <  t  <  T  (11) 


By  the  Wiener-Khintchine  theorem,  we  know  that 

R(T)  "J/"S(W)  elTW<iw 

S(w)  ■  f  r(t)  e"iTWdT 


(12) 


S(w)  is  the  power  spectral  density  of  the  random  sound  pressure 
p(t).  We  now  suppose  that  S(«)  is  a  rational  function  of  w  so 
that  it  can  he  written  as  the  ratio  of  two  polynomials  in  w. 


s(.)  -  iilislfl 

D  ((iw)2) 


(15) 


D  is  at  least  one  degree  greater  than  I.  The  argument  of  the 
N  and  D  polynomials  was  chosen  as  a  square  of  w  because  of  the 
general  requirements  placed  on  S(«)  as  a  power  spectral  density 
function.  We  wish  to  show  that  when  the  process  has  such  a 
spectral  density  as  (13),  the  functions  t  (t)  satisfy  a  linear 
constant  coefficient  differential  equation. 


We  start  hy  inserting  (13)  into  (3)  to  obtain 
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K(m)2)  n 


(t*  )dwdt' 


l\l2V« 


Let 


K(i»)2) 


p  >  q. 


(HO 


(15) 


a  and  bfc  are  constants, 
with  the  operator 


Sow  operate  upon  both  sides  of  (Ik) 


kmO 


2k 

d _ 

*k  2k 
dt 


to  give 


1.  rV 

2 irjy  J 


dumV^'H 

K(1«)  )  “ 

■  l>ln|2  t>td?)  *n(t) 


(t' )dt*dw 


(16) 


Multiplying  by  H((i»)2)  under  the  Integral  J 
operating  on  the  integral  with  “* 
written  as 


is  equivalent  to 
Thus  (l6)  can  be 


~&f0  [,{3P)/e'  “*->  Vt,)dt' 

2  (17) 

-  I\l* 

But 

rm  i(t-v)w 

J  e  dw  -  2t r  B(t-t* )  (18) 
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vhere  6(t-t’ )  is  the  Dirac  delta  function.  Using  (l8)  in  (17), 
ve  find 

*n(V)dt'  ■  |XJ2  “‘a?1  *n(t) 
or 


fhia  is  the  basic  differential  equation  that  must  be  satisfied 

by  the  expansion  functions  t  (t). 

n 

The  solution  of  (19)  will  contain  A  and  p  arbitrary 
constants  A  .  When  this  solution  is  substituted  back  in  the 
integral  equation  (ll),  it  will  be  seen  that  the  integral  equation 
cannot  be  satisfied  except  for  certain  values  of  the 
characteristic  values,  and  the  constants  must  satisfy  certain 
conditions.  These  conditions  show  that  for  each  n  there  is  only 
r»n^  independent  constant.  This  one  constant  is  to  be  determined 
from  the  normalizing  conditions  on  +  (t),  namely. 


/ 


tn(t)  tj  (t)dt  -  {  1  n  “  J 


0  n  ^  j 


(20) 


The  procedure  Just  outlined  will  now  be  carried  out  explicitly. 
The  differential  equation  (19)  becomes 


q/2  ,2k 


I 


k«0 


.  ,/r  d2\(t) 

**  at2* 


k  2k 
dt 


■  i\i’i 


k-0 


(21) 


Try  a  solution 


.  x  .  at 
t_(t)  »  A  e  n 
n  n 


(22) 


A  is  a  constant  and  a  is  a  parameter.  Putting  (22)  in  (21) 
ve  find  that  0  must  bs  a  root  of 
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(2?) 


q/2  p/2 

I  vf  ■  kj*i  v: 

k-0  k-0 


2k 

n 


Bils  equation  will  have,  in  general  (special  situations  can  be 
taken  care  of),  p  roots  £  Q  -  ♦  0  .  (X  ),  k  »  1>..,p/2.  Half 
of  the  roots  are  simply  the  negatives  or  the  other  half.  A 
general  solution  of  (21)  is  then 

P/2 

*n(t)  -  X  '“■**  +  V  (2M 

kul 

The  autocorrelation  function  is  given  by 


R(t) 


i.  r’m  i")2) 

21  i  Kd-)2) 


iru 

e  du 


(25) 


This  integral  can  be  evaluated  by  contour  integration  using 
semicircular  arcs  and  yields 

P/2 

f  y  v*“k’  T  t‘t’>0 

k-1 


are  the  roots  of 

D((iw)2)  *  0  (2Yj 


lying  in  the  upper  half  of  the  complex  plane  as  shown  in 
Figure  2 
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There  are  p  roots  of  (27)  +  wk>  k  *  !>•••■>  p/2» 

B,  are  given  by 
k 


B 

k 


N((i«k)2) 


dP <(i“)2) 


dw 


U) 


UL 


(28) 


Inserting  (2k)  and  (26)  into  (11),  ve  obtain 
p/2  P/2 

r’y  -*«*,«.•  . 


rCv"-1-- ■!  <.  v 


-  -  -Onkt’ 


)  dt’ 


P/2 


I* 


+  at  -  -a ,  t, 

enk  +  A  ,e  nk  ) 
n'  Z_u.  nk  nk 

k« 
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Carrying  oat  the  integrations  we  obtain 


yyBA*  £il_  .  yyBA- 

V  4-  *  **  “nk'1".  v  4  *  *  '“ok'1 


Q  -i 4* 
nk  s 


-JSv 


nk  -Q  +lw 
nk  s 


_yy  BA*  pi.  _yyBA- 

V  4  •  ■*  “nk’1".  V  4-  *  “*  -“nk'1". 


i?v 

•Hv 


+  e(flnk+i‘*r)T  _iw  t 


nk  a  ,  +  iw 
nk  s 


-  «(_ank  ♦  i"*)1  t 


nk  -0  .  ♦  iw 
nk  s 


e  s 


l\'2J‘<k-  "nk 


<Mx  +  a^  •'“■>**) 


This  equation  oan  be  satisfied  identically  in  t  only  if  the 
sum  of  the  coefficients  of  the  different  functions  of  time 
are  equated  to  zero.  This  gives  the  system 


(a) 


0>) 


P/2  A+  A- 

I  Vs-^r  ♦  1  *•-** 

k51  nk  8  nk  s 


P/2  A+  e(Qnk  + 


!*• 


nk 


a  .  ♦  iw 

nk  s 


A  e 
nk 


(-Onk  + 


-0  .  +lw 
nk  s 


) 


P/2 


<c)  I A*  -  H7^r}  B. 

Sal  nk  s  nk  s 


1> • • •  ,p/2 

|x  |V. 

1  n1  nk. 


k  ■  1, .  •  • >  p/2 


(29) 
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P/2 


L 


Av{“ 
nk  Q 


nk  s 


-a 


nk 


1  ■  }  B  ■  |X  |  A  k  ■  1, ...,p/2 

•flu  s  1  n'  nk. 


Equations  (29a)  and  (29b)  are  a  system  of  p  linear  homogeneous 
equations  for  the  determination  of  the  A^"  and  A  •  For  a  non¬ 
trivial  solution,  the  determinant  of  the  coefficients  must  vanish, 
namely, 


-1 


-1 


0  ,-iu,  0  ,-iu 

nl  1  nl 


a  #  +  iw 

ni  1 


°ni+  ^1 


-1 


a  ,  -iu 

nl  2 


a  .  +  iw 

ni  1 


1 

V  -%/2 


e 


(a  +  iu  )t 

no  o 


a  .  +  iw 

nl  1 


(30) 

-  o 


(a  iu  )t 

e  nl  o 

a  ,  +  iw. 
nl  1 


(-Ono  ♦ 

O  _ 

-a  .  ♦  iw 

ni  1 


From  this  equation  ve  determine  the  eigenvalues  A^,  n  ■  1,2, . . . 
This  requires  that  ve  know  the  functions 


(31) 


These  functions  f  (A  )  are  determined  by  solving  for  the  roots 
of  (23)  in  terms  of  ft  .  We  can  do  this  analytically  in  terms 
of  algebraic  functions  for  equations  of  fourth  or  lower  degree. 
This  would  mean  that  problems  could  in  principle  be  handled  for 
D((iw)2)  Qf  eighth  or  lower  degrees. 
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Equations  (29c)  and  (29d)  will  be  found  to  be  identically 
satisfied.  All  constants  A^v,  k  »  1,  ..  p/2-1,  A^,  k  -  1,..., 
p/2  can  be  solved  for  In  terns  of  Aj^/2  using  equations  (29a)  and 
(29b).  A^p/2  i*  determined  from  the  normalizing  condition 


(52) 


The  case  where  the  power  spectral  density  Is  a  rational 
function  of  frequency  is  a  very  important  one  and  covers  a 
wide  range  of  practical  cases. 

b)  Orthogonal  Expansion  of  Processes  with  Markoff 
Autocorrelation  Functions 


Suppose  that  the  autocorrelation  function  of  a  stationary 
random  acoustic  noise  process  p(t)  of  mean  zero  is  given  by 


E  (p(t)  p(t  +t)}  -  R(t)  -  A  e“^T^ 

2 

A  -  E(p  (t)}  ■  R(o),  p  »  constant 


(33) 


This  function  is  shown  in  Figure  3» 


-pItI 

Figure  3.  Graph  of  the  Autocorrelation  Function  Ae^ 
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The  spectral  density  corresponding  to  (33)  1a 


S(«) 


e-eM  -  im 


dr 


2Ap 

2  2 
-(**)•* 


function  Is  plotted  In  Figure  4. 


(34) 


Figure  4.  Spectral  Density  Function  of  a  Process  with  Markoff 
Autocorrelation  Function 

From  (34)t  ve  see  that  the  process  under  consideration  has  a 
rational  spectral  density  function  vith 


H((lu)2)  -  apA 
d«i»)2)  -  -(M2  +  e2 

The  differential  equation  (19)  becomes  in  this  case 

2fSA*  (t)  -  |A  |2[-  P2*  ] 

“  ”  dt2 

Let 


u 


2 

n 


-P2  + 


2gA 


TSien  (3 6)  can  be  written  as 


(35) 


(36) 

(37) 


l4 


(38) 


d2t 


n  2 

T  +  «  ♦  -  0 

2  n  n 


If 


dt 

P2<  141011  (3®)  has  solution 


lu  t 


.  i.i  jjui  (,  -in  t 

♦  (t)  ■  a  e  n  +  b  e  n  ,  u  >  0. 
n  n  n  n 


(39) 


an  bn  •**  artltnu7  constants.  Substituting  (33)  into 
the  Integral  equation  (11),  we  find 


X  I 

A/  •*\(t,)dt»  +aJ  ePTtn(t')df  -  l\|2*n(t),  t  -  t-t» 


Using  (39)  and  carrying  out  the  Integrations  gives 


isfet  -pt  .-Hit  -pt 

A*  — - - -  ♦  At  5 - - -  + 

n  _  .  n 

P  +  iw  p  -  iw 

n  n 

P(t-T)  +  HT  P(t-T)  -  luuT  -  -iu_t 

♦  Aa  - — -  -  +  Ab  — - - 

n  ^  +  lu  n  ^  .  in 

n  n 


l\l*  ♦  |x_|2 1..-^ 


n'  n 


This  equation  is  identically  satisfied  in  t  if 


Aa 


Aa 


p  +  iw  -e  +  iw 


-  •V  an 


Ab 


Ab 


P  -  i«  -p  -  iw 


|A  |2  b 
n1  n 


(*0) 
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-Aa 

n 


Ab 

n 


P  +  iw 


P  -  iw 


(-6  +  iw  )T 
Aa  ev  n' 

n _ 

-p  +  iw 


Ab  +  *“n>T 

_ n _ 

P  +  iW 


0 


Bie  first  two  conditions  in  (to))  are  identically  satisfied.  She 
second  two  are  'tec  linear  homogeneous  equations  for  a  and  b  . 
For  a  nontrivial  solution,  we  require  n 


1 


1 


P  +  iw 
n 


(-P  +  iw  )T 
e  n 

-p  + 1** 


-(P  +  iw  )T 
"6  n 

P  +  iw 
n 


-  o  (ki) 


This  reduces  to 


w  cos  ^  w  T  +  p  sin  ^  w  T  »  0 
n  2  n  2  n 


Equation  (to>)  determines  the  eigenvalues  A  .  It  is  clearly 
a  complicated  equation  and  must  be  solved  numerically. 

From  (to)) 


b 

n 


p  -  iw 
_ n  a 

p  +  iw  n 
n 


(*3) 


Therefore 


l£ 


♦  (t) .  .  (.‘“i* .  e-%t 

n  n  p+l«n 


] 


0+i«  '  n 

n 


(« cob  w  t  +  p  sin  w  t} 

u  n  ^ 


<W) 


From  the  normal  lying  condition  we  find 


1‘nl2  ■ 


p  +  lu 


f  •{’ 


-1 


I «  cos  «  t  ♦  p  sin  w  t|  dt)  (45) 


Ofte  orthogonal  expansion  of  p(t)  now  takes  the  fora 

f  l\l*n 

p(t)  ~L  £T7S  (“n  co*  V  +  P  ““  VJ  Xn 


(Mi) 


n»l 


Tb*  xn  be  determined  when  a  realized  value  of  p(t) 

is  given.  This  in  principle  completes  the  problem. 

c)  Other  Methods  for  Solving  the  Integral  Equatl on 

Another  method  of  solving  the  Integral  equation  (11)  much 
used  by  physicists  is  to  expand  the  kernel  R(t-t')  into  a  suit¬ 
able  set  of  orthoraoimal  functions.  By  this  means  the  problem 
c*°  be  reduced  to  the  problem  of  finding  the  non-zero  roots  of 
an  infinite  determinant.  By  Judicious  choice  of  the  expansion 
this  process  can  be  made  tractable.  (Ref.  5.) 

d)  Other  Solutions 

1.  Stationary  Band  Limited  Gaussian  White  Noise  of 
Mean  Zero  (Ref.  4.) 

The  expansion  interval  in  -.  <  t  <  «  .  The  autocorrelation 
function  for  the  process  is 


R(t)  -  <p2> 


sin  2ttBt 

2T7BT 


and  is  shown  in  Figure  5. 


t  -  t-t* 


m 
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Figure  5.  Autocorrelation  Function  of  Band  Limited  White  Noise 

B  *  u^~  U)1  18  the  handvldth  of  the  noise.  The  spectral  density 

appears  as  in  Figure  6. 


Figure  6.  Spectral  Density  of  Band  Limited  White  Noise 


Note 

P(t) 


that  \p  >»  S  B.  The  proper  functions  for  the  expansion  of 
on  the  interval  (-•,  •)  are  found  to  be 


♦n(t)  -  45 

2irB(t  -  n/2B) 


m 
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)  la  an  orthonormal  set  but  t  not  solutions  of  Eq,.  11. 

n  n 

The  expansion  takes  the  form 

xt>  -  to®  -1  f  **Msa^*m  ,  m 

xk-m  271®  (t  -  n/zB) 

where 

i  -as  tRft)'1  [lw  ^  ™LvV?»), .  dt  <*» 

2uB(t-V2B) 


2.  Band  Limited  Stationary  formal  White  lolse  on  a 
finite  Internal.  (Ref .6.) 


Ihe  correlation  function  is  still  of  the  form  (Vf) ,  but  the 
expansion  Is  desired  on  the  interval  (o,T).  The  proper  functions 
and  eigenvalues  are 


*»<*> 


8*1*  [l*T,  2/T(t-l)] 
on 

Vi- 

on 


o  <  t  <  * 


X  -  f  p2>  T(rJ^(BT,  l)]2  (51) 

n  on 

•  (*aB>2 

*on  "  2X  2J  4-  1 
J 

R  (c ,y)  and  3  Ac,  cos  9)  are  respectively,  the  prolate 
nl  ni 

spheroidal  Bessel  function  and  prolate  spheroidal  lagendxe 
function  as  given  by  Norse  and  Fesbbech  (Ref.  5) .  These 
functions  have  been  tabulated.  The  expansion  reads 


P(t)  -^VT  tn(t)  Xn,  o  <  t  <  T 
Xn  *  Jr/ P(t)  ♦n(t)4t 


(52) 
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zero 


X  are  Independent  nornal  variates  of  mean 


J.  White  Molse 

Correlation  and  Spectral  density  functions  are 


R(t)  -  Sq6(t) 

S(«)  -  so 


(55) 


Sq  is  a  constant.  The  Integral  equation  in  this  case  becomes 


T 

S  6(t-t' )  f  (t’Jdt* 
o  n 


I\|“  *,W 


0r  2 

S  ♦  (t)  -  |x  |  ♦  (t) 

on  n  n 


This  can  be  satisfied  by  any  orthonormal  set  of  functions  if 

we  take  X  ■  */s  for  all  n.  The  expansion  becomes 
no 


p(t)  «/s  )  t  (t)  X 

oZ-.  n  n 

T 

Xn  ■  jr  /  *<*>  V‘)dt 


(54) 


The  results  for  white  noise  are  very  Interesting  and 
potentially  very  useful.  We  can  make  the  general  statement 
that  very  broad  band  random  signals  can  always  be  approximately 
expanded  into  orthogonal  series  using  any  suitable  set  of 
orthonormal  functions. 


IV.  Methods  of  Approximating  Acoustic  Hoise  Fields 
a)  Optimal  Approximation  Techniques 

Once  X  and  t  (t)  are  determined  from  solving  the  integral 
equation,  $  can  §e  obtained  when  a  typical  recording  p(t)  Is 
available.  °We  get  a  set  of  realized  values  for  the  X  .  The 
series  )  X  tn(t)  ■  p(t)  Is  then  an  analytic  representation  of 

the  realized  function  p(t).  The  sirens  could  be  used  to  try 
to  simulate  this  series  or  the  typical  recording. 
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In  general,  when  we  hare  an  analytic  expression  for  p(t),  we 
need  a  synthesis  procedure  to  approxlaate  p(t)  with  the  sirens. 
We  sight  try  to  approxlaate  p(t)  itself  In  the  tine  domain  by 
requiring  the  neaa  square  difference  between  the  sirens'  output 

Ipk(t)  and  p(t)  to  be  a  minimum.  This  requirement  Is  achieved 
by  minimizing 


|  / ?p(t)  -^pk(t»2dt 


(55) 


with  respect  to  the  adjustable  parameters  of  the  sirens  such  as 
center  frequencies,  amplitude  of  outputs,  percent  modulation  etc. 
Let  be  a  parameter  of  the  kth  siren,  then  this  synthesis 

procedure  requires 


I 

■I  J,k<t»2dt 

k»l 


0 


k 


1;2 p  e • e  p  H 


*k 


1,2, ..., 


(56) 


IT  Is  the  number  of  sirens  and  L  Is  the  number  of  adjustable 
parameters  of  the  k^1  siren.  ($6)  leads  to  a  system  of  equations 
for  the  optimizing  parameters.  Unfortunately,  the  system  Is 
nonlinear  In  the  parameters  In  almost  all  cases  and  It  Is  difficult 
to  solve  even  with  digital  computers.  Here  we  try  to  optimally 
synthesize  one  realized  value  of  p(t).  Since  p(t)  Is  a  sample 
function  from  a  random  process,  it  may  be  better  to  approximate 
some  statistical  characteristics  of  the  process. 


Other  synthesis  criteria  might  be: 


1.  Minimize  mean  square  difference  between  power  spectral 
density  of  p(t)  and  that  of  X’k<'>- 


2. 

3. 


Optimize  probability  density  of  p(t). 
T  N 

Minimize  |  E {[  [p(t)  -£pk(t)]2dt} 
®  k-1 


with  respect  to  the  parameters  of  the  sirens.  We  assume  here  that 
some  of  the  parameters  of  sirens  are  random  in  character. 


b)  Spectral  Broadening  of  a  Single-Tone  Siren. 
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Spectral  broadening  of  a  single-tone  siren  can  be  achieved 
by  various  fonts  of  modulation  or  by  switching  techniques.  Let 
us  consider  first  spectral  broadening  by  modulation. 

1.  Amplitude  Modulation  (Deterministic) 

A  single-tone  siren  output  may  be  represented  by 


p(t)  -  A  cos  («ct  +  f)  (57) 

where  A  is  the  amplitude  of  the  output,  w  is  the  center  fre¬ 
quency  of  the  siren,  and  f  is  a  phase  angle.  Such  a  siren  output 
as  (57)  contains  only  one  spectral  component.  Suppose  by  some 
means  we  are  able  to  vary  A  as  a  function  of  time.  In  particular, 
let 

A  •  A(t)  -  A  (1  ♦  A  cos  w  t)  (58) 

o  m 

A  ,  A,  and  u  are  constants 

°  4  UT’U/1‘ 

den  (57)  reads  ^  Met 

CM)  CeJ)  tect  CVJ  7  ~ 

P(t)  -  A  (1  +  A  cos  »  t)  cos  (at  +  f)  (59) 

o  me 

or 

p(t)  »  A  cos(w  t+f)+^AA  cos  [(«  +  w  )t  +  •]  + 
o  c  2  o  cm 

1  (60) 

+  —  A  A  cos  [(w  -  id  )  t  +  •] 

2  o  '  c  m  T 

By  this  modulation,  we  see  that  p(  t)  now  has  three  spectral 
components.  It  is  easy  to  see  that  if 

U 

p(t)  ■*„! (i  +  \  cos  w^t)  cos  («ct  +  f)  (61) 

k-1 

p(t)  will  have  2H  +  1  spectral  components.  Equation  (6l)  includes 
a  vide  variety  of  amplitude  modulation  such  as  arbitrary  periodic 
modulation  and  almost  periodic  modulation.  We  note,  however,  that 
as  long  as  A(t)  is  a  deterministic  function  of  time  the  spectrum 
is  discrete. 
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2.  frequency  Modulation  (Deierninistic) 

The  instantaneous  angular  frequency  is  defined  a a 

“<*>  *  fj  <“/+  ♦(*))  («) 

irtiere  the  phase  angle  f  ia  now  considered  a  function  of  tine. 

If  im  wary  the  frequency  by  a  single  tone,  the  siren  output 
takes  the  font 

p(t)  .  Ao  eo.[«ct  +  A-  .la  (63) 

m 

A  and  art  constants 

It  is  well  known  that  this  can  be  written  in  the  fora 

m 

P(t)  .  |  Ao  co.  act  +  |*J  (1  ♦  (.!)»)  )  (#0.  (U'%  +  ^)t 

,  (6A) 

♦  cos  (w  -  nu  )t) 

e  a 

fWa  which  we  see  that  there  are  an  infinite  nunber  of  discrete 
spectral  components.  is  the  Bessel  function  of  the  first 

kind  of  order  n.  B*  general  appearance  of  the  intensity  spectra 
of  p(t)  for  snail,  interne dlate ,  and  large  values  of  A  /«  are 
shown  in  figure  7  with  their  envelopes.  * 


«5 


Tor  arbitrary  periodic  or  almost  periodic  frequency  modulation, 
the  siren  output  is  given  by 


p(t) 


A  cos 
o 


.  y  \ 

{«  t  +  )  - 

c  L>  « . 


sin  w  t) 

mk 


k-1 


m 


a  u  are  constants. 

Tc  mk 

•The  number  of  sidebands  increases  tremendously  as  R  Increases, 
but  the  intensity  spectral  is  still  discrete.  Only  for  an 
aperiodic  save  can  ve  obtain  a  continuous  intensity 

spectrum  or  spectral  density. 


3.)  Amplitude  Modulation  ( Random) 

Continuous  spectral  broadening  of  a  single-tone  siren  can  be 
achieved  through  random  modulation.  Consider 

p(t)  -  A  [1  +  A  P  ft)]  cos  w  t  (66) 

O  A  c 

vhere  P  (t)  is  a  stationary  normal  random  process  and  A^,  A,  and 

u  are  constants.  13*  spectral  density  of  p(t)  for  a  modulating 
noise  with  spectral  density  function  given  by 

s  (.)  .  s  (67) 

M  o 

S  and  u  »  constants 
o  b 


is  shown  in  Figure  8  for  various  modulation  indices. 

It  is  seen  that  spectral  boradening  can  be  quite  considerable. 
Rote  that  finite  energy  is  in  the  carrier  for  all  modulation  Indices, 
shown. 


4-.)  Frequency  Modulation  (Random) 

In  this  case  we  have 

p(t)  ■  A^  cos  twct  +  P^(t)] 

where  P^(t)  is  a  stationary  normal  random  process.  Aq  and  «c  are 

constant.  Bie  spectrum  of  p(t)  for  ?M(t)  taken  to  have  spectral 
density 


(68) 


t A 


V-> 


* 

1  ♦  m2/B2 


if  <P„>  6 

p '»  (  8„  +  "*'*  )  «9) 

3n7 


i»  shown  In  Figure  9  for  rariou*  conditions.  B  is  the  band- 

vldth  of  the  noise.  <TP  Si«  V  , 

\rM  Z1*  fflean  square  value  of  nodulating 


B? 


[se^J 


isoo 
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Figure  8.  Spectral  Density  of  Noise  Modulated 
Single -Tone  Siren 
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A  tO  h 


Figure  9.  Spectral  Density  of  a  Single-fcne  Siren  FM  by  Stationary 
Horaal  Boise  With  Markoff  Autocorrelation  Function 
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If  |i  <  <  l,  the  spectral  density  is  very  narrow  about  «c 
(little  broadening) .  For  very  large  values  of  n,  the 
spectral  density  is  bread  about  «e. 

5.)  Spectral  Broadening  by  Random  Switching 

Consider  a  stationary  noise  pressure  of  the  type  shown  in 
Figure  10. 


Figure  10.  Baadoa  Square  Have 


The  zero  crossings  occur  at  purely  random  instants  t^,t  ,.... 
Such  an  output  would  be  produced  (except  for  a  carrier  term 
cos  w  t)  by  an  ideal  siren  that  vaS  being  switched  on  and  off 
at  purely  random  times.  The  spectral  density  function  of  the 
pressure  in  Figure  10  is 


8  (a)  -  “ 

N  i r 


1  +  (-5  ) 

when  tj^-tf  is  a  Poisson  distributed  random  variable. 

If  we  use  I(t)  as  an  amplitude  modulation  function  of  a  single- 
tone  siren  output,  we  have  the  sound  pressure 


(70) 


P(t)  -  H(t)  cos  [wt  +  f] 


(71) 


^  is  a  random  phase  uniformly  distributed  between  o  and  2f. 
autocorrelation  function  of  p(t)  is 


88 


Rp  -  E{M(t)  cos  [*ct  +  f]  f  H(t  +  t)  cos  [«c(t  +  T)  +  f])  (72) 

i  r2r 

Rp  •  E(H(t)  W(t  +  t)}  ~ J  cos  [»ct  +  f)  cos  [w  t)  +?  )  df 

o  c 

Rp(T)  •  |  (t)  COS  fc»cT 


spectral  density  is  obtained  by  use  of  the  Wldner-Klntchlne 
theorem.  We  find 

sp(«)  Rp(t)  e”WTit 


V“> .  .-<•  «•.*$  VT)dt 

SpM  "  £  (V— c}  +  SE(W  +  Wc)} 


(73) 
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V.  Conclusions  and  Recommendations  for  Further  Study 


The  energy  of  a  » ingle -tone  siren  can  be  spread  over  a  fairly 
large  frequency  bend  about  the  center  frequency  of  the  siren  by 
various  sorts  of  modulation.  The  most  effective  way  of  spreading 
the  sound  energy  over  the  frequency  band  is  by  random  amplitude 
or  frequency  modulation • 

A  given  acoustic  noise  field  may  be  approximated  (in  spectral 
density)  by  centering  the  sirens*  center  frequency  at  those  which 
constitute  the  given  signal  or  by  dividing  the  spectral  density 
of  the  given  signal  into  equal  energy  bands  and  assigning  an  energy 
•t>and  to  each  siren.  The  approximation  can  be  improved  by  amplitude 
•nfl  frequency  modulating  each  siren.  Approximation  of  a  given 
acoustic  noise  field  in  some  optimal  way  leads  to  the  solution 
of  a  system  of  non-linear  equations. 

Further  work  may  be  carried  out  to  find  other  solutions  to 
the  integral  equation  occurring  in  the  Karhunen-Loeve  expansion 
of  the  random  noise  field.  It  would  be  particularly  desirable 
to  find  some  exact  solution  so  that  the  proper  expansions  could 
be  compared  with  a  Fourier  expansion  of  the  field.  Another 
direction  in  which  a  continued  investigation  might  bear  fruit  is 
that  of taying  "to  devise  methods  of  solving  the  nonlinear  equations 
arising  in  trying  to  optimally  approximate  some  characteristic  of 
the  given  acoustic  noise  field  with  the  sirens.  Still  another 
area  worth  looking  into  further  is  spectral  broadening  by  randomi¬ 
zation  of  the  sirens'  parameters. 
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